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Massless Dirac fermions in graphene at charge neutrality form a strongly interacting system in 
which charge-neutral modes comprised of correlated particle-hole excitations play a dominant role. 
The neutral modes are essentially decoupled from charge dynamics in the absence of a magnetic field, 
but become strongly coupled to charge modes when a field is applied. We show that these ideas 
explain the recently observed giant magnetodrag, arising in classically weak fields when electron 
density in the layers is tuned near charge neutrality. We predict strong Hall drag in this regime, 
which is in stark departure from the weak coupling regime, where theory predicts the absence of 
Hall drag. Hall drag and magnetodrag arise in a wide temperature range and at classically weak 
magnetic fields, and feature an unusually strong dependence on field and carrier density. 



Graphene near charge neutrality (CN) hosts an in- 
triguing charge-compensated electron-hole system with 
unique properties. Theory predicts a number of inter- 
esting phenomena for this system: strong interactions 
and criticalitypH!], momentum-conserving collisions giv- 
ing rise to finite resistivity in the absence of disorder j5j- 
7 ; strong thermopower resulting from coupling be- 
tween charge modes and neutral modes [5111 lj; anoma- 
lous thermodynamics [3] and new collective modes [T2"lll3| . 
Our understanding of the behavior at CN would benefit 
tremendously from introducing ways to efficiently couple 
the novel neutral modes predicted at CN to charge modes 
which can be easily probed in transport measurements. 
There is a long history, particularly in areas such as 
plasma physics and astrophysics, of employing magnetic 
field for such a purpose. The hydrodynamics of charge- 
neutral plasmas is ultra-sensitive to the presence of exter- 
nal magnetic fields. Further, dynamically generated mag- 
netic fields can result in striking new phenomena, such 
as thermomagnetic waves, predicted in astrophysics fB] 
and later discovered in a 3D semimetal|15j. 

A new promising system in which magnetotransport 
at CN can be probed are atomically thin graphene- 
based heterostructures [TrJl IT?] . In these systems, strong 
Coulomb coupling between adjacent layers leads to 
strong interlayer correlations. A striking manifestation 
of these correlations is provided by Coulomb dragfPTl- 
125] . in which current applied in one (active) layer in- 
duces a voltage in the adjacent (passive) layer. Strong 
interactions manifest themselves in a drag resistance that 
peaks near CN [T7] . Interestingly, the measurements [T7] 
revealed dramatic magnetic field dependence, with the 
drag resistivity at CN increasing by more than an order 
of magnitude (and changing sign) upon application of a 
relatively weak B field. In particular, for fields as low as 
B ~ 0.1 T the longitudinal drag at CN features a five- 
fold increase as compared to drag at B = [TT]- Strong 
magnetic field dependence of drag has been observed pre- 
viously in other double layer two-dimensional electron 
gas (2DEG) heterostructures [26H28] . however these ex- 
periments were carried out in the quantum Hall regime, 
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FIG. 1: Magnetodrag in a double layer graphene heterostruc- 
ture in the presence of a finite magnetic field, (a) Schematic of 
charge current, temperature gradients, and electric field in the 
two layers. (b,c) Predicted den sity dependence of magneto- 
drag resistance, i?jj rag , from Eq.(12| and Eq.(|l4|. Parameter 

values: B = 0.6 T, n = 10 11 cm^% T = 150 K, and p = ^ . 
The B = dependence taken from the model of drag at zero 
B field [241 125] . (d) Experimentally measured magnetodrag re- 
sistivity in G/hBN/G heterostructures, reproduced from Ref. 
[17] for the same B values as in (c). Application of magnetic 
field leads to a giant increase of drag at CN. Note the simi- 
larity between data and theoretically predicted drag density 
dependence, B dependence, and sign. 



whereas the anomalous magnetodrag found in Ref. [17] 
occurs at classically weak fields B < 1 T. This puzzling 
behavior, which is not yet understood, will be explained 
below (see Figjlj). 

Another outstanding question that magnetodrag mea- 
surements may help settle is that of Hall drag. It has long 
been debated whether a Hall voltage can arise in the pas- 
sive layer in the presence of a DC current in the active 
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layer [2^1 [3U] . In particular, theory predicts a vanishing 
Hall drag [29 , 30 . This is because transferred momentum 
is parallel to velocity allowing only a longitudinal "back- 
current" to develop in the passive layer. As discussed 
below, this conclusion is valid only in the weak coupling 
regime. At strong coupling, charge and energy transport 
become dominated by long-range neutral modes, which 
in our system is described by electron-lattice temperature 
imbalance. At CN, the magnitude of the cross-couplings 
between charge and energy currents becomes large, pro- 
ducing a finite Hall drag 

v H = R d ™ s i l{ . (i) 

Our theory predicts i?^ ag values which are consistent 
with observations |31j but orders of magnitude greater 
than those observed in earlier experiments (see below). 

The flow of charge and energy, which are completely 
decoupled in the absence of magnetic field (in a spatially 
uniform system), become strongly coupled in the pres- 
ence of a magnetic field. As we will see, the effects aris- 
ing due to such coupling dominate Coulomb drag at CN. 
Heat current and an electric field, induced by charge cur- 
rent and temperature gradients, are coupled via the ther- 
moelectric effect altered by the B field, 

VT 

j q = Qj, E = Q — . (2) 

Here Q is a 2 x 2 matrix, of which diagonal compo- 
nents describe the Pelteir and Thompson effects, and off- 
diagonal components describe the Nernst-Ettingshausen 
effect. Crucially, vertical energy transfer between layers 
does not discriminate between longitudinal and trans- 
verse heat currents since temperature profile is a scalar 
field. This stands in contrast to conventional momentum 
driven drag, where momentum transfer is parallel to the 
applied current [2^1130] ■ As illustrated in Fig. [T^,, when a 
longitudinal charge current is applied in the active layer 
(for B ^ 0) a transverse (Ettingshausen) heat current 
develops in both layers through efficient vertical energy 
transfer. Nernst voltage in the passive layer results in 
a longitudinal magnetodrag (see Fig. [T|3,c). The result- 
ing density and field dependence strongly resembles the 
observed magneto drag 1 1 7] . reproduced in Fig. [TJi. 

As we will see below, this mechanism also leads to Hall 
drag. As illustrated in Fig(T|and Fig[2| we find Hall and 
magnetodrag resistances, -R^ ag and thai are large 

and peak near CN. As we will see, these large values 
arise even for classically weak fields B ~ 0.1 T. Indeed, 
both i?^ ag and i?j| lag exceed by two orders of magni- 
tude those found previously in GaAs/AlGaAs 2DEG het- 
erostructures [27Tff2"5| at similar fields. Additionally, both 
quantities exhibit sharp growth as a function of B. The 
steep differential magnetodrag resistance can be used in 
sensitive field and temperature sensing applications. 

In our analysis we will focus on a Hall-bar geometry 
of two parallel rectangular layers of dimensions L x W, 



pictured in Fig[T^. For the sake of simplicity, we treat 
the electric and heat currents as independent of the x 
coordinate along the bar. In layer 1, current is injected 
at x = —L/2 and drained at x = L/2. In layer 2, the 
Hall drag voltage arising across the device, Vh and the 
longitudinal drag voltage, Vy , are evaluated as 




E™dy (3) 



The electric and thermal variables may depend on the 
transverse coordinate y, reflecting thermal coupling be- 
tween the two layers. 

To obtain the electric field in layer 2 induced by cur- 
rent applied in layer 1, we first need to understand the 
coupling of temperature profiles T 12 (r) in the two layers. 
Energy transport in the system can be described by 

-VkjWTi + a{5T x -5T 2 ) + A<57\ = -V • (Q (1) j) 
-Vk 2 V5T 2 + a(8T 2 - 5T X ) + XST 2 = (4) 

where a is the energy transfer rate between the two layers 
|25j . A is the electron-lattice cooling rate, and STi — Ti — 
Tq. Here To is the lattice temperature of the two layers. 

The boundary conditions are given by the condition of 
electric current being tangential to the free boundaries 
at y — ±W/2, and zero temperature imbalance at the 
Hall-bar ends, x — ±L/2, since the current and voltage 
contacts act as ideal heat sinks. The electric current 
parallel to the boundaries y = ±W/2 gives rise to the 
Ettingshausen heat current that may have a component 
transverse to the Hall bar. The divergence of this heat 
current, appearing on the right hand side of Eq.Q, acts 
as an effective boundary delta function source in the heat 
transport equations. 

We consider the case of a spatially uniform Q in both 
layers. The ideal heat sinks at x — ±L/2 mean that no 
temperature imbalance develops in the x-direction (ex- 
cept for some "fringing" heat currents near the contacts 
which give a contribution small in W/L, we will ignore 
these fringing currents in the following discussion). In- 
deed far away from these contacts, no temperature gra- 
dients are sustained in the x-direction and we reduce 
our problem Eq. Q to a quasi-one dimensional prob- 
lem with temperature profiles that only depend on the 
y-coordinate for W < L. As a result, the only heat source 
arises from the Ettingshausen effect Q^j = (Qyxj)y- 

To describe transport in the presence of such a source, 
we will expand temperature variables in both layers in 
a suitable orthonormal set of functions. Here it will be 
convenient to use eigenstates of the operator az with zero 
Neumann boundary conditions at y = ±W/2, given by 

My) = Am (^ y ), „„(„)= A si n(^±Hy), 

A = (2/W) 1 / 2 , n = 0, 1,2... From the symmetry of the 
source in Eq.Q we expect STi^iy) to be odd in y. Thus 



3 



only the functions v n (y) are relevant, giving 
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For each n we obtain a pair of algebraic equations 

qlmSTi + a{5T x - ST 2 ) + XS^ = F n 
qlK 2 Sf 2 + a{5f 2 -Sfi) + XST 2 = 



(5) 



where k% j2 



.(1,2) 



and F n = 2A(-l) n Q$j. Solving 



Eq.([5]), we find the temperature profile in layer 2: 

aF n 



ST 2 (y) 



^ Li(q n )L 2 (q n ) - a 2 



v n {y), 



(6) 



where Li{q n ) — Kiq 2 + a + X (i = 1,2). Since electron- 
lattice cooling is very slow [33J [33] , with the correspond- 
ing cooling length values in excess of few microns, we 
will suppress A in what follows. Because the boundaries 
in the transverse (y-direction) are free, a finite temper- 
ature imbalance between the edges can arise, given by 
AT = ST 2 (y = W/2) - ST 2 (y = -W/2). We find 



(1) i 



AT = 4A 2 aL £yx] _ o \ ~» ^syx 

^ LlL2 -a 2 wk^qia+eaiy 



(7) 



where we defined k — k\ + k 2 and a length scale £ = 
\J K\K 2 /ak. We evaluate the sum using the identity 
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Connecting AT with the drag voltage, and in particu- 
lar determining its sign, requires taking full account of 
Onsager reciprocity. This analysis is presented below. 

In the same way that the applied charge current, j, 
in layer 1 causes a heat current (Pcltcir/Ettingshausen), 
a temperature imbalance in layer 2, AT, can sustain 
voltage drops across the sample (Thermopower/Nernst). 
These two effects are related by Onsager reciprocity con- 
straints. The cross couplings in the coupled energy and 
charge transport equations |34j arise from 
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L22 , 




(9) 



where L are 2x2 matrices and e is the carrier charge. 
In this notation, the electrical conductivity equals g_ — 
e 2 L 11 /T, and thermal conductivity is k — L 22 /T 2 . Com- 
paring to the heat current due to an applied charge cur- 
rent, Eq.pl), we identify L 21 = — eQL n . 
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FIG. 2: Hall drag in a graphene heterostructure in the pres- 
ence of a magnetic field, (a) Schematic of charge current, 
temperature gradients, and electric field in the two layers. 
(b,c) Density dependence of Hall drag resistance, predicted 
from Eq. ( 12 1 and Eq. ( 14 1 for the same parameter values as in 



Fig[T] (d) Density dependence of Q X x, Qxy, see text. 

Onsager reciprocity demands that the cross-couplings 
obey L 12 (-B) = L 2 r 1 (— B) where B is the applied magnetic 
field (note the transposed matrix). In an isotropic sys- 
tem the off-diagonal components of L obey jj xv ^ (£?) = 
Lr yx \—B). As a result, Onsager reciprocity reduces to 



L 12 (B)=L 21 (B) 



(10) 



in an isotropic system. Applying Eq. [I0]to Eq. [9] in an 
open circuit, we find E = -e _1 V// = T^L^QLuVT. 
For an isotropic system Q = Q xx l — iQ xy <r 2 , L = L xx l — 
iL xy a 2 , so that [Q,L] = 0, which gives Eq.|2j). 
Using Eq.((3]) and Eq.Q we obtain 

(11) 

giving the magnetodrag and Hall drag resistance values 
-G(0, 




R 



drag 
H 



Tk 



-0 (1) (2) 

^i X y ^Cx X ■> 



R 



drag 



- LG (Z) n(i)n(2) 
WTk ^v^y 



(12) 

Here we have used the fact that in an isotropic system 
Q xx = Qyy and Q xy — —Q yx . When £/W <C 1, we have 
G — >• 1 so that saturates to become independent 

of the the energy transfer rate between the two layers, a. 
However, when £/W ^> 1, we have G — > so that R^^ 
vanish with decreasing energy transfer rate, a. 

To describe the density and B field dependence, we use 
a simple model for Q. Measurements indicate [HI IH] that 
thermopower and the Nernst effect in graphene are well 
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described by the Mott formula |3S], giving 



Wms — „ / ,Pmk~ 
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(13) 



where p is the magneto-resistance and u is the chemi- 
cal potential. Here we concentrate on classically weak 
magnetic fields relevant for recent experiments |17) . We 

model the tensor p — ^ P p H as follows: 



P\\ = 



Po 



-Bn 



n 2 /rig 



Ph 



e(n 2 + no) 2 ) 



(14) 



where po is the resistivity peak at the Dirac point, n 
is the carrier density, and parameter uq accounts for 
broadening of the Dirac point due to disorder. We ac- 
count for disorder broadening of the density of states, 



J/ 2 



dn/dp = {n 2 + n 2 ) 1 ' i {2/{k% 2 v 2 f )) 

Using Eqs.Q,0,|l4j) we obtain i?f| rag (in Fig{TJj,c) 

and i?^ ag (in Fig|2p^c). In that, we used the parameter 
values no = 10 11 cm -2 , po = 3^2 , and a representative 
temperature, T = 150 K. These values match the char- 
acteristics (disorder broadening, hq, and peak resistivity, 
po) of G/hBN heterostructures described in Ref.|17j. As 
a sanity check, we plot the components of Q (in Figj2ji) 
which show the behavior near CN matching thermopower 
and Nernst effects measured in graphene [5J [S]. 

Plotting magnetodrag in Fig. jlj), we find i?jj lag peaks 
at dual CN showing a giant and negative trans-resistance. 
At its peak, Ru lag is two orders of magnitude larger than 
magnetodrag measured in GaAs quantum wells [2l)H2"5] 
where drag reisistivities of several ohms were seen at sim- 
ilarly weak fields. Magnetodrag in Fig. [TJ; exhibits a 
steep B dependence, -Rjj r pf ak oc — B 2 , that bears a strik- 
ing resemblance to measurements in Ref.|17j reproduced 
in Fig[ljl. In particular, our model explains the negative 
sign of the measured magnetodrag peak. 

Hall drag, shown in Figj2]D,c, is large and sign- 
changing. Finite i?^ ag is a 'smoking gun' of strong 
coupling to the neutral modes at CN in graphene (see 
Fig. [2]}). Already for classically weak magnetic fields, 
B < 0.1 T, this strong coupling manifests in a giant i?^ ag 
of several hundred ohms (Fig. [2J2). This is two orders of 
magnitude larger than the few-ohm Hall drag resistance 
measured in GaAs quantum wells [25] at similar fields. 
i?^J ag also exhibits dependence on magnetic field, with 
peak heights increasing linearly with B (Fig(2j;). Inter- 
estingly, only carrier density in layer 2 controls the sign 
of i?^ ag , breaking the n\ •<->• 712 symmetry between layers 
- a unique signature of Hall drag in graphene. 

In summary, magnetic field has a strong impact on 
drag at CN because it induces strong coupling between 
neutral and charge modes, which are completely decou- 
pled in the absence of an applied magnetic field in a 
uniform system. Field-induced mode coupling leads to 



giant drag that dwarfs the conventional momentum drag 
contribution as well as a remnant drag due to spatial 
inhomogeneity[25]. Our estimates indicate that these two 
contributions are orders of magnitude smaller than the 
predicted magnetodrag, which also has an opposite sign. 
The giant magnetodrag and Hall drag values attained at 
classically weak magnetic fields, along with the unique 
density dependence and sign, make these effects easy to 
identify in experiment. The predicted magnetodrag is 
in good agreement with findings in Ref.jTTj. Magnetic 
field, coupled with drag measurements at CN, provides a 
unique tool for probing the neutral modes in graphene. 
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